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We observe an obvious anomalous line-shape of the e'*'e~ -^ hadrons total cross sections in the 
energy region between 3.700 and 3.872 GeV from the data samples taken with the BES-H detector 
at the BEPC Collider. Re-analysis of the data shows that it is inconsistent with the explanation for 
only one simple ^(3770) resonance with a statistical significance of 7a. The anomalous line-shape 
may be explained by two possible enhancements of the inclusive hadron production near the center- 
of-mass energies of 3.764 GeV and 3.779 GeV, indicating that either there is likely a new structure 
in addition to the ^(3770) resonance around 3.773 GeV, or there are some physics effects reflecting 
the DD production dynamics. 
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In the energy range from 3.700 to 3.872 GeV, the well 
established ^'(3770) resonance is believed to be the only 
observed structure. This resonance has been identified 
to be a mixture of Z?-wave and S'-wave of angular mo- 
mentum eigenstates of the cc system. In addition, the 
■0(3770) resonance is expected to decay into DD me- 
son pairs with a branching fraction that is greater than 
98%. However, there is a long-standing puzzle in the ex- 
isting measurements of '(/'(3770_) production and decays. 
Before recent BES-II 0, H S i, i] and CLEO-c re- 
sults published, existing data indicated that about 38% 
of ■!/'(3770) does not decay to DD final states [7|. Re- 
cently, the BES Collaboration measured the branching 
fraction of ■!/'(3770) decays to DD to be 5[-0(377O) -^ 
DD] = (85 ±5)% d S, and directly measured 
the non-DD branching fraction of ■(/'(3770) decayto be 
B[4>{3770) ^non-Dp] = (13.4 ± 5.0 ± 3.6)% [3 and 
B[^l3770) ^non-DD] = (15.1 ± 5.6 ± 1.8)% under 
assumption that there is only one simple i/'(3770) res- 
onance in the energy region between 3.700 and 3.872 
GeV. In the last two years, the BES and CLEO Col- 
laborations have searched for exclusive non-DD decays 
of -0(3770). However, the summed non-DD branching 
fraction measured by both the BES and CLEO Collabo- 
rations remains to be less than 2% [l|, la| . To understand 
why the measured inclusive non-DD branching fraction 
is substantially larger than 2%, in addition to continu- 
ing searching for more possible non-DD decay modes of 
■0(3770), it is worth going back to carefully examine the 
previous measurements of the ■0(3770) parameters. 

An examination of analyses previously reported by the 
BES Collaboration in Refs. [2,l9| shows that the fits to the 
observed hadronic cross sections or R values are rather 
poor for the fine-grained energy scan cross section mea- 
surements [see Fig. 4(a) in Ref. [3| and Fig. 1 in Ref. [9|] 
even though the branching fraction for 0(3770) ^non- 
DD was left as a free parameter in the fits. In this letter, 
we present a reanalysis of the observed inclusive hadronic 
cross sections to better understand the hadronic annihi- 
lation structure in the energy region between 3.700 and 
3.872 GeV. 

The measurements of the observed inclusive hadronic 
cross sections are discussed in detail in the Refs. [2, 3, 9, 
|lO|. The observed inclusive hadronic cross sections ob- 
tained from the cross section scan data taken in March 
2003 and in December 2003 are illustrated in Fig. [1] [17] 
by dot with error bars, where the error bars are the 
combined statistical and point-to-point systematic un- 
certainties. The systematic uncertainty includes the sta- 
tistical uncertainty of the luminosity, the uncertainties of 
the Monte Carlo efficiencies for detections of the Bhabha 
scattering events and the hadronic events, as well as the 
uncertainty of the observed cross sections due to the re- 
producibility (±0.1 MeV) of setting the BEPC machine 
energy. The cm. (center-of-mass) energy of the BEPC 
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FIG. 1: The measured inclusive hadronic cross sections versus 
the cm. energy for the two data sets taken in March and 
December 2003; the fit is done with two incoherent amplitudes 
(solution 1), see text for detail. 



0(3686) and J/^p- The measured masses of 0(3686) and 
J/ip at BEPC are obtained by analyzing 6 data sets of 
0(3686) scan and 2 data sets of J/tp scan performed dur- 
ing the time periods of collecting the finer cross section 
scan data. The uncertainty in the calibrated energy for 
the combined two finer cross section scan data sets to- 
gether is about ±0.5 MeV. 

A close examination of the energy region (from 3.74 
to 3.80 GeV) around 3.777 GeV shows that the slopes 
of the observed cross sections on the two sides of the 
peak are quite different; with the slope of the high energy 
side of the peak substantially larger than that of the low 
energy side. It conflicts with the expectations for only 
one resonance in this energy region, since the effects of 
the initial state radiation (ISR) and the DD production 
threshold as well as the energy dependence of the DD 
scattering amplitudes due to the Blatt-Weisskopt bar- 
rier [llj would all make the slope at the high energy 
side of the peak less steep relative to the slope on the 
low-energy side. This anomalous shape seen in the pre- 
cision measurement indicates that one simple resonance 
hypothesis is quite questionable to fit the current data. 
Instead of the conventional definition of the ^0(3770) de- 
cay width T{Ecm), if the dynamics of DD scattering or 
some reasonable model describing the DD scattering can 
give some special form of F (E'en!) and mass shift for which 
the scattering amplitude gets zero or node at the rather 
low D meson momentum (Pp ~ 0.4 GeV) to adapt the 
unusal decline around 3.8 GeV in the cross section line 
shape, the anomalous line-shape of the cross sections for 
e+e^ -^ hadrons might be understood. 

However, as shown in this work, it can not be excluded 



to the "0(3770) resonance in the energy region between 
3.700 and 3.872 GeV, which and its interference with the 
■0(3770) ampHtude distort the Une-shape of the observed 
cross section from that expected if there was only one 
resonance in the region. 

To investigate whether there are some new structures 
in addition to the -0(3770) resonance in the energy region 
between 3.700 and 3.872 GeV, we fit the observed cross 
sections with one or two amphtudes in the energy region. 
The expected cross section cr^^^'^^ (Ecm) consisting of four 
components can be given as 



with 
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■1 • -1 expect / T-i \ expect/ 771 \ expect / j-i \ 

m which Cr^4377o)(£^cm), CT^/^ (EJem), ^.0(3686) (^"n)' 

and cr]^™{Ecin) are, respectively, the expected cross 
sections for i?s(3770) -^ hadrons, J/0 — + hadrons, 
0(3686) —^ hadrons, and continuum light hadron pro- 
duction at the cm. energy Ecm, and i?s(3770) denotes 
the full structure around 3.773 GeV. The expected cross 
sections are obtained from the Born order cross sections 
for these processes and the ISR corrections [12|, [l3| . 

For the i?s(3770) resonance(s), we use one or two 
pure P-wave Breit-Wigner amplitude(s) with energy- 
dependent total widths [3, S to fit the observed 
hadronic cross sections. The two amplitudes are expected 
as 



^j(-C'cm) — 



i27rr|'=rhad 



{El^~Mf)+trf'iEcm)M, 



(J = l,2), 



(2) 

where Mj, T^^, T^^'^, and r^°*(s) are the masses, leptonic 
widths, hadronic widths, and the total widths of the two 
resonances, respectively. r*°*(£^cm) is chosen to be en- 
ergy dependent [3, S |9|. For two amplitude hypothesis, 
concerning the possible interference between the two am- 
plitudes, we use two extreme schemes to see if we can get 
better description for the anomalous line shape. In the 
first scheme, we ignore the possible interference; and in 
the second, we assume the complete interference between 
the two amplitudes. These two schemes give the Solution 

1 and Solution 2, respectively. The Born order cross sec- 
tion for i?s (3770) production in Solution 1 and Solution 

2 can, respectively, be written as 



and 



^fl.(3770)(^cm) - |Ai(Ecm)P + |A2(Scm)|' (3) 



^fls(3770)(^cm) = |^l(i?cm) + e'^As (i?c„0 P , (4) 



where the is the relative phase difference between the 
two amplitudes. 

The non-resonant background shape is taken as 
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where o'Lt™(i?cm) is the observed cross section for light 
hadronic event production given in Refs. [2, Isl, cr^^ _ (s) 



is the Born cross section for e^e^ —^ fi'^ii~ , pjjo aiidpD+ 
{E^o and £'£>+) are the momenta (energies) of Z?° and 
D~^ mesons produced at the nominal energy ^/s, 9oo and 

6*-! are the step functions to account for the thresholds 

of the D^D^ and D^D" meson pair production, respec- 
tively; / is a parameter to be fitted. The effect of energy 
spread on the observed cross sections is also considered 
in the analysis. 

In the following, ignoring the tiny difference of the de- 
tection efficiencies determined from the different schemes 
as described above, we fit the observed cross sections pre- 
sented in Fig. [1] and in Fig. [2l respectively, with the ex- 
pected cross sections given in Eq. (1) in two schemes. In 
the first case, it is defined in Eq. (3) and the fits give 
the results of the Solution 1. In the second case, it is 
defined in the Eq. (4) and the fit gives the results of the 
Solution 2. As a comparison we also fit the cross sections 
with the conventional one Briet-Wigner form of 0(3770) 
resonance as the definition of the i?s(3770) for the one 
resonance hypothesis. In the fits, we fix r = 1.5 fm (r 
is the interaction radius of the cc system) [2|, y, |9[ and 
fix the J/-0 parameters at the values given in PDG07 [8|]; 
the -0(3686) and i?s (3770) resonance parameters are left 
free, i?uds and / @, Q are also left free. 

As shown in Fig. [T] and in Fig. [21 the circles with error 
bars show the observed cross sections. The red lines in 
both of the figures and in the sub-figures (a) inserted in 
Fig. [T] and Fig. [2] represent the fitted values of the cross 
sections of Solution 1 and Solution 2. The green lines 
in the sub-figures (a) show the fit to the observed cross 
sections for the one amplitude hypothesis. The circles 
with error bars in red as shown in the sub-figuares (b) 
inserted in Fig.[l]and Fig.[2]show the measured net cross 
sections, which are obtained by subtracting the contri- 
butions from J/tp and -0(3686) decays to hadrons, the 
continuum hadron production and the interference term 
of the two amplitudes in _Rs(3770) definition for the Solu- 
tion 2; the blue lines show the fit to the net cross sections 
from the two resonances for both of the Solution 1 and 
Solution 2, respectively. 

The 2nd, the 3rd and the ith columns of Table H] sum- 
marize, respectively, the results of the fits for the Solution 
1 and the Solution 2 of the two amplitude hypothesis, and 
for the one amplitude hypothesis, where the first errors 
are from the fit and the second systematic. For the mea- 
sured masses, the second errors mainly arise from the 
uncertainty of the BEPC machine energy calibration for 
the combined two data sets together. For the one res- 



TABLE I: The fitted results for the data taken in March 2003 and December 2003. 



Quantity 


two amplitudes 


two amplitudes 


one amplitude i/)(3770) and G(3900) amplitudes 




(without interference) 


(interference) 


(interference) 




Solution 1 


Solution 2 


Solution 3 



X /ndof 



M 



[MeV] 

r!,^,«««, [kev] 



i/>(3686) 



- 1^(3686) 
-'ee 
V(3686) 



[keV] 



Ml [MeV] 
rl°' [MeV] 
rr [eV] 
Mz [MeV] 
r'2°' [MeV] 

or aa [MeV] 
TT [eV] 

or C 

f 



125/103 = 1.21 
3685.5 ±0.0 ±0.5 

312±34±1 

2.24 ±0.04 ±0.11 

3765.0 ±2.4 ±0.5 

28.5 ±4.6 ±0.1 

155 ± 34 ± 8 

3777.0 ±0.6 ±0.5 

12.3 ±2.4 ±0.1 

93 ± 26 ± 9 



0.4 ±5.6 ±0.6 



112/102 = 1.10 182/106 = 1.72 

3685.5 ±0.0 ±0.5 3685.5 ± 0.0 ± 0.5 

311±38±1 304±36±1 

2.23 ±0.04 ±0.11 2.24 ±0.04 ±0.11 

3762.6 ± 11.8 ± 0.5 3773.3 ± 0.5 ± 0.5 

49.9 ±32.1 ±0.1 28.2 ±2.1 ±0.1 

186 ± 201 ± 8 260 ± 21 ± 8 
3781.0 ±1.3 ±0.5 
19.3 ±3.1 ±0.1 

243 ± 160 ± 9 



(158 ± 334 ± 5) 
5.2 ±2.5 ±0.6 



0.0 ±0.5 ±0.6 



170/104 = 1.63 

3685.5 ±0.0 ±0.5 

293 ± 36 ± 1 

2.23 ±0.04 ±0.11 

3774.4 ± 0.5 ± 0.5 

28.6 ±2.3 ±0.1 

264 ± 23 ± 8 

3943.0 (fixed) 

54 (fixed) 

0.243 (fixed) 
(150 ±23 ±5) 
0.0 ±1.2 ±0.6 



parameters as listed in the Ath column of Table [H These 
measured values of the resonance parameters are consis- 
tent within error with the world averages 14 1 [l8| and 
with the earlier BES measurements P [S'l obtained by 
analyzing the two data samples separately. The fit gives 
the mass difference between the ■4>{3770) and ip{3686) 
resonances to be Am = 87.8 ± 0.5 MeV. However, the 
large y^ /ndof in the Ath column of Table |T] gives the 
fit probability of less than 7 x 10^^, meaning that the 
one resonance hypothesis is strongly incomparable with 
the present precision measurement data. On the con- 
trary, the x^ change for the two hypotheses in Solution 
1 is (182 — 125) = 57 with a reduction of 3 degrees of 
freedom. This indicates that the signal significance for 
the new structure is 7.0(7. The x^ change for the two 
hypotheses in the Solution 2 is 70 with a reduction of 
4 degrees of freedom. This indicates that the statisti- 
cal significance of the new structure is 7.6(7. Comparing 
the fits for the Solution 1 and Solution 2, we find that 
the x^ change of 13 with a reduction of 1 degree of free- 
dom. The significance of the interference between the two 
Breit-Wigner amplitudes is 3.6(7, which indicates that the 
two amplitudes likely interfere somehow with each other. 
The actual situation of the interference would be some- 
where between the two cases. It depends on what are the 
exact final states of the possible new structure decays. 

However, it is noted that the fitted value / = 5.2±2.5± 
0.6 in the Solution 2 would lead to a huge DD production 
cross section at higher energy region and there exists an 
evident dip of the inclusive hadronic cross section around 
-Ecm=3.80 GeV. These indicate that, instead of only the 
continuum DD production, there might be a broad struc- 
ture whose peak is at higher energy than 3.83 GeV and 




3.7 3.75 3.8 3.85 

^cm [GeV] 

FIG. 2: The observed inclusive hadronic cross sections versus 
the nominal cm. energies for the combined data sets taken in 
March and December 2003; the fit was done with two coherent 
amplitudes for 7?s(3770) (Solution 2). 



BELLE collaborations [16,] observed 0(3900). To con- 
sider the effect of the G(3900) on the observed cross 
sections, instead of the first two solutions for the two 
structure hypotheses one may adopt the third approach 
by including the new component of DD production am- 
plitude of G(3900). The fitting procedure is analogous 
to Solution 2. However, the amplitude A2{Ecm) in Eq. 
(4) is replaced by a square root product of a parame- 
ter C and a Gaussian function G. The mass and the 
standard deviation of G are fixed at the measured val- 



ing to the DD cross section as the one measured by 
BABAR at 3.943 GeV. The red hne in Fig. [^a) rep- 
resents the fitted values of the cross sections, which is 
obtained from the fit under assumption that the ■0(3770) 
and G(3900) amplitudes interfere with each other; the fit- 
ted value from the hypotheses for only ■!/'(3770) amphtude 
(blue line), from Solution 1 (yellow line) and from Solu- 
tion 2 (green line) are also illustrated in Fig. [3ja). The 
5th column of Table U summarizes the results (Solution 
3) of the fit including G(3900). The fit gives a rather 
poor fit probability of less than 5 x 10~^, which does 
not significantly improve the fit from the one resonance 
hypothesis. If we consider three coherent amplitudes 
in the fit by replacing \Ai{Ecm) + e^'^A2(£'cm)P with 
\AiiE,m)+e'*'A2iEem)+e'^'G{E,^)\^ inEq. (4), where 
G is the G(3900) structure, we obtain almost the same 
results as these shown in Solution 2 in Table I instead of 
/ = 5.2 ± 2.5 ± 0.6. This fit gives / = 2.7 ± 6.4 ± 0.6, 
which is comparable with the inclusive hadronic cross sec- 
tion measurements at the higher energy region. Fig.[3l^b) 
shows the ratio of the residual between the observed cross 
section and the fitted value for the one 0(3770) ampli- 
tude hypothesis to the error of the observed cross section. 
The variation of the ratio with Ecm indicates that there 
is more likely some new structure additional to 0(3770) 
resonance. 
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FIG. 3: (a) the observed inclusive hadronic cross section 
versus the nominal cm. energy; (b) ratio of residual to error 
of observed cross section; (see text). 

In summary, by re-analyzing the line-shape of the cross 
sections for e^e^ — > hadrons, we find that it does not 
describe the cross section shape well with the hypothese 
that only one simple ■0(3770) resonance exists in the en- 
ergy region from 3.700 to 3.872 GeV. If there are no other 
dynamics effects which distort the pure D-wave Breit- 
Weigner shape of the cross sections, the analysis shows 



one simple 0^(3770) resonance there at 7a statistical sig- 
nificance, indicating that there might be evidence for a 
new structure additional to the single ^0(3770) resonance. 
However, if there are some dynamics effects distorting the 
pure D-wave Breit-Weigner shape of the cross sections, 
such as the rescattering of DD leading to the significant 
energy dependence of the wave function in the DD de- 
cays of the 0(3770) resonance, one has to consider those 
effects in the measurements of the resonance parameters 
of ^0(3770), since these effects would definitely shift the 
measured values of the resonance parameters. Anyway, 
the large non-DD branching fraction of 01(3770) decays 
measured previously [2, l3| may partially be due to the 
assumption that there is only one simple resonance in 
the energy region between 3.700 and 3.872 GeV in the 
previous measurements of the ^"(3770) parameters. 
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